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In this paper, we present a general formula for obtaining 
the reduced density opeator for any biparticle pure entangled 
state. Using this formula, we derive, in a compact form, the 
explicit formula of the entanglement for any bipartical pure 
entangled Gaussian state. In the case of Gaussian states, the 
criteria of separabelity can be naturely obtained by the for- 
mula. For non-Gaussian states, we also show the usefulness 
of the method presented in this paper. 

PACS: 03.65.Ud, 42.50.-p 



Quantum entanglement is one of the essential features 
of quantum mechanics. It has been an interesting and im- 
portant topic since it was first noted by Einstein, Podol- 
sky, and Rosen (EPR) [Q. Now, it has attracted much at- 
tention because of immense progresses both in the foun- 
dation of quantum mechanics and in the burgeoning field 
of quantum information theory. Issues such as the rela- 
tionship between entanglement and quantum nonlocality 
and the violations of Bell's inequality for both discrete 
|||,D and continuous variable states || further enrich the 
contents of quantum mechanics. In recent years, quan- 
tum entanglement is also viewed as a useful "resource" for 
various kinds of quantum information processing. The 
original EPR states of continuous variables can be cre- 
ated either by the nondegenerate optical parametric am- 
plifier (^-0], or simply by using passive optical elements 
(e.g., beam splitters) |8|-|l0|. While most quantum infor- 
mation protocols were initially developed for quantum 
systems with discrete quantum variables, quantum in- 
formation processing based on quantum states with con- 
tinuous variable has also been proposed, e.g., quantum 
teleportation pd| , quantum error correction p2| , quan- 
tum computation |e| , entanglement purification p4 | and 
cloning [ p"5| |. Thus, inseparability criteria jl^-[T^] and en- 
tanglement quantification |l9| ] in continuous variable sys- 
tems has become an issue of practical importance. 

The measurement of entanglement for any biparticle 
pure entangled state, i.e., the von Neumann entropy of 
either partial trace of the density operator for the state, 
is considered as a "good measurement" of quantum en- 
tanglement. Parker et al. |19| have developed an elegant 
method to calculate the entanglement of the biparticle 
pure entangled states with continuous variables by means 



of the integral eigenvalue equations in the coordinate- 
momentum space. In fact, their method is based on the 
Schmidt decompositions. But in certain cases, it is some- 
what harder to find such a "Schmidt basis" for the contin- 
uous variable systems. In this paper, we propose an alter- 
native method to calculate the entanglement entropy for 
continuous variable states in Fock space. With the help of 
linear quantum transformation theory (LQTT) pO| , we 
obtain an explicit formula of the entanglement entropy 
for any Gaussian state. With this formula, we can eas- 
ily calculate the Gaussian states' entanglement without 
knowing the eigenvalues of the reduced density operator 
for the system. Furthermore, from the derivation of the 
formula, we naturally obtain the necessary and sufficient 
condition of separability for a biparticle pure entangled 
state with continuous variables. For some non-Gaussian 
states, one can also calculate the eigenvalues of the re- 
duced density operators, from which the entanglement of 
this bipartite system may be easily obtained. 

It is well known that the entanglement entropy of a 
biparticle pure entangled state is defined by 



tr px In pi — - tr p 2 In p 2 , 



(1) 



where the reduced density operator pi(2) is a partial trace 
of density operator pi 2 for the state, i.e., 



Pl(2) 



tr 



2(1) Pl2- 
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In order to give a general formula of reduced density op- 
erator pi for any biparticle system, we use the following 
formula MM 



n 



: exp 



(a\, a. 



+ (9 Zl ,d Z2 ) 



x (Z\Q\Z) 



(3) 



where £1 is an arbitrary bosonic operator, a] (i = 1,2) 
is the usual bosonic creation operator for the subsystem 
i; \Z) = \Z\) ® \Z 2 ) denote the usual coherent states, 
: : means the normal ordering and z — means ^1(2) — 

Z l(2) 

p\ 2 for the biparticle system as 



we write the density operator 



P12 



exp 



(a\,al) 



a. 



{d Zl , d z 



(Z\ Pl2 \Z) \ z=0 



(4) 



1 



where \Z) = \Z\,Z%). With the help of overcomplete 
of coherent states J S-^z- \Z 2 ) 22 (Z 2 \ = 1, and taking the 
partial trace over subspace of the subsystem 2, one can 
get the reduced density operator p\ as follows 



Pi = tr 2 P12 = tr 2 
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(at,4)( /D+ia^a^d) 



x : e 
d 2 z 
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J ~~ : eX P ( a l^i* + 



x exp 2 2 *5 Z . + d Z2 z 2 * (Z \p l2 \ Z) 



(5) 



In the following we will consider two classes of states: 
the Gaussian states and the non-Gaussian ones. Let us 
first consider the following Gaussian states 



p 12 = A : exp < i 



(oUi)Mi(° x 



+ (4, a 2 )M 2 f ° 2> ) + 2(4, oi)Mia f° 2 

\ a 2/ \02 



(6) 



where An is a normalization factor, (i = 1,2) is Her- 
mitian matrices, and 



M12 = 



e / 
/* e* 



with e and / being two arbitrary complex numbers. In 
Eq. (|^) the linear terms of a and arc not included 
since they do not affect the entanglement of the biparticle 
Gaussian states |lJJ] . In this case, the matrix elements of 
pi 2 with respect to \Z) is 



{z* 1 ,zx)M 1 [ ^) +(zlz 2 ) 



xM 2 A +2(^,zi)M 1: 

\z 2 J \z 2 

Using the following Gaussian integration formula 

d 2 z ( 1 fz*\ ( z* 

— expj--(z*,z)Q( \ +(u,v)[ 

[- detQl^expj^u.^Q^r 



(7) 



(8) 



where Q = Q is nonsingular and u (v) is an arbitrary 
complex number, and the formula 

exp(z 2 *<9 z j + <9 Z2 z 2 ) *(zi,zi;z 2 ,z 2 ) \ Z2=Z * =0 
= *(z^,zi;z2*,Z2) , (9) 
we have, after substituting Eq. (Q) into Eq. (||), 



An 



Pi = 



V- det Af 2 ' 
According to LQTT 



KaI,a 1 )(M 1 -M 12 M- 1 M 12 )(^l) 



M = 



if one denotes 
a <i 



(10) 



(11) 



there is a map of Af defined as 

D(M) = (Mi - Mi 2 M 2 ~ 1 Afi 2 )E i3 1 
„-i i E, 



cr 1 - 1 C -M ^ 



1 - C" 



fl2) 



where complex constants a, 6, c and d are determined by 
the known matrices Mu2) an( i M12 in Eq. (|iTj|). Then 
we can immediately rewrite the normally ordered form 
in Eq. © as @ 



pi = /I exp 



oi 



(13) 



where negative Hermitian matrix iV = InAf, J^b 

1 A ■, A ■ 

, and A is constant 



-1 



A = A 



- det M 2 

Substituting Eq. ((TJ) into Eq. (Q), we get 



in which 



-A 



Z(P)\nA-—Z(0) 



0=-i 



Z(/3) = tr exp | - 1 [a\ , a x )NY> B 
Using the result in Ref. ]22] ] 

Z(/3)=trexp|-^(4,ai)iVE B ^| 
= |det(e^-l)r^ , 



(14) 



(15) 



(16) 



direct calculation yields 



d/3 



Z(P) 1/3= 



1 



dct(e 



-N 



1) 2 tr- 



N 



1 - e 



jv • 



(17) 



(18) 



Substituting Eq. (|T^) into Eq. (15) we get the explicit 
expression of the entanglement for the state in Eq. (ph 



E = —A |det(e /3Ar - 1 



, , 1 N 

lnA+ 2 tr T^ 



(19) 



Thus, from Eq. ( p9| ) one can conveniently obtain the 
entanglement of states in Eq. M) only by calculating the 



2 



eigenvalues of negative Hermitian matrix N, which can 
be easily obtained. This is the main result of this paper. 
It is worth pointing out that from Eq. ([l9]), the necessary 
and sufficient condition of separability for the states in 
Eq. (H) naturally reads 



In A = - tr ■ 

2 e 



N 



N 



1 



(20) 



As an application of Eq. jl9| ) , we consider the entan- 
gled states produced by a beam splitter 0] 



M 13 = B(M)Si(Ci)Sa(Ca)|00>, 
where the beam splitter operator B(9 : 4>) is J23| 

B(0,< 



exp 



(a\a 2 e tv - a^e - ^) 



and the single-mode squeezed operator S(C) is |24 



5(C) = exp 



i(C« 2 -C« t2 ) 



(21) 



(22) 



(23) 



Making use of LQTT, we can write the normally ordered 
form of the density operator in Eq. (Ell) as follows 



P12 



)i2i2 0>\ = M : exp | i 



+ (4,a 2 )A/ 2 f a ^ +2( a t 1 , ai )M 12 f^ 



(24) 



where 
Mi = 
Mi 2 = 



a 1 

1 a* 

5 
S* 



Mo = - 



(3 1 
1 /?* 



a = i^t tanh | Ci I cos 2 + e 2lip ^- tanh |C 2 1 sin 2 I 

ICil IC2I 



P = e" 2 ^^ tanh |(i| sin 2 + tanh | Ca | cos 2 0, 

ICil IC2 1 

5 = 1 sin (20) (A tanh |C 2 | - ^ tanh | Ci I e"*). 

2 K2I ICil 



<2 



Using Eqs. © and (0) yields 
1 

I <5 1 cosh |Ci cosh| IC2I 



.4 



M 



a <i 
6 c 



(25) 



with 



0= i(l + 6d), c = - — 0- 



d = 



L (<5 2 /3* + a(l - |/3| 2 )) , 6 



= -d* 



Finally, we get the entanglement of |'0)as m ^1- (^1) 
1 



# = 



\S\ cosh I (1 cosh| |^2 1 
x [ln(|J|cosh|Cicosh| |C 2 |) 

1 + A 



■In A 



(27) 



2(1 -A) 

where A is any of the roots in the following equation 

A 2 - (a + c)A + 1 = 0. (28) 

We can consider a special case of Eq. (^lj) by choosing 
= f, Ci = s t e^ (i = E2) and <p = l f (I = 0,1,2 ) 
0. In this case Eq. @ is just the result of Ref. 0, 



£ = cosh 2 \s\ ln(cosh 2 \s\) - sinh 2 \s\ ln(sinh 2 |s|). (29) 

where s = \{s\e tip + s 2 e~ l¥ '). Furthermore, if we choose 

(p = 0, C2 = — Ci = r an( l = f j we get Ai = tanh 2 r. we 
then get the well known result for the two-mode squeezed 
vacuum state, whose entanglement is [[19| 



i? = cosh 2 r ln(cosh 2 r) — sinh 2 r ln(sinh 2 r) 



(30) 



For non-Gaussian states, generally speaking, their en- 
tanglement can not be written in the compact form as 
in Eq. (U9). But for some states such as SU(2) coher- 
ent states(see, e.g., 0]) and the photon-added coherent 
states p5[ , the entanglement can be easily obtained by 
solving the eigenvalue problem of the reduced operator of 
Eq. (|3|). To show this, we consider the following example 

Pi 



= B\n 1 ,n 2 ) = e e ^ a2 - a ^\n 1 ,n 2 ) 



(31) 



The density operator p\ 2 of this system can be written, 
by following the above procedure, as 



P12 = B\ni,n 2 ) (n 2 ,ni|-B t 
f 1 



1 d ni d n2 
m\n 2 \ da ni d/3™ 2 



x : exp 



aia 1 )M 1 ( a i)+(ala 2 )Mj a2 

\ a lJ \ a 2 



+2(4,ai)M 



12 



a 2 



\a=/3=0 



(32) 



Here 



Mi = (a cos 2 9 + (3 sin 2 6 - l) 01, 
M 2 = (/3 cos 2 + a sin 2 - l) a u 
M3 = sin cos 0(/3 — a)(Ti , 

(26) where o\ is Pauli matrix. Using Eq. (^|), we derive 



(33) 



3 



d" 2 



PA 



n v n 2 \ da n * df3 n i (^cos 2 9 + asm 2 9 - l) 



a[3 — (a cos 9 + (3 sin 2 6») 



\ (f3cos 2 9 + asm 2 9- l) J 

Thus we get the eigenvalues as follows 
1 d ni d" 2 1 



|q=/3=0 



(34) 



A 



Ni 



n v n 2 l da n ± d(3 n ^ (3 cos 2 (9 + a sin 2 9 

-iJVi 



1 



a/3-1 



/? cos 2 # + a sin 



|a=/3=-l (35) 



Note that in Eq. ( j35| ) we take the replacements a = a +1 
and P = (3 + 1 . Finally, we get the entanglement of the 
state in Eq. Q 



A'. 



In A 



AT,- 



(36) 



In summary, we have presented a general formula of re- 
duced density operator for any biparticle pure entangled 
state of continuous variables. Based on this formula, we 
have derived the explicit expression of the entanglement 
for any Gaussian pure state in a compact form. From the 
examples mentioned above, it is easy to see that Gaus- 
sian states' entanglement can be conveniently calculated 
by the formula given in this paper. For the case of non- 
Gaussian states, we have also shown the usefulness of the 
formula for calculating the reduced density operator, as 
given above. In this sense, the method to calculate the 
entanglement for biparticle pure entangled state is easier 
to operate. 
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